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Abstract 

Measurements continuous in time were consistently introduced in 
quantum mechanics and applications worked out, mainly in quantum 
òptics. In this context a quantum filtering theory has been developed 
giving the reduced state after the measurement when a certain trajec- 
tory of the measured observables is registered (the a posteriori states). 
In this paper a new derivation of filtering equations is presented, in 
the cases of counting processes and of measurement processes of dif- 
fusive type. It is also shown that the equation for the a posteriori 
dynamics in the diffusive case can be obtained, by a suitable limit, 
from that one in the counting case. Moreover, the paper is intended 
to clarify the meaning of the various concepts involved and to discuss 
the connections among them. As an illustration of the theory, simple 
models are worked out. 
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Introduction 



Usually in quantum mechanics only instantaneous measurements are con- 
sidered, but by using the notion of instrument [H| also measurements 
continuous in time were consistently introduced E| [E] and applications 
worked out El EE21 US] HH1 ; see also [201 ESI - 

Now a natural question is: if during a continuous measurement a cer- 
tain trajectory of the measured observable is registered, what is the state of 
the system soon after, conditioned upon this information (the "a posteriori" 
state)? By using ideas from the classical filtering theory for stochastic pro- 
cesses and the formulation of continuous measurements in terms of quantum 
stochastic differential equations [131 EU UH 1221 121] , a stochastic equation for 
the a posteriori states has been obtained [21] |2E]· The main purpose of this 
paper is indeed that of clarifying the meaning of that equation by present- 
ing a natural derivation of it in terms of instruments, independently from 
any notion related to quantum stochastic calculus, and by discussing some 
models. Moreover, we shall discuss the connections among various things 
appeared in the literature about what can be called a quantum version of 
the theory of stochastic processes (continuous measurements) and filtering 
theory (a posteriori states). 

Let us start by recalling the important notions of instrument and of a 
posteriori states. The notion of instrument has been introduced in the op- 
erational approach to quantum mechanics pp. Let a quantum system be 
described in a separable Hilbert space TC and denote by B(7í) and T(7i) the 
Banach spaces of the bounded operators on 7í and the trace-class operators, 
respectively. Let (í), E) be a measurable space (Çl a set and E a cx-algebra of 
subsets of Q). An instrument |I] [H] 2T is a map from E into the space of the 
linear bounded operators on T(7í) such that (i) Z(B) is completely positive 
[2H| for any BeE, (ii) J2j%(Bj)g = X( Ui Bj)g for any sequence of pairwise 
disjoint elements of E and any g in T(7í) (convergence in trace norm), (iii) 
Tr{l(Q.) 6 } = r L·{ e }, \/g G T(H). 

The instrument X is an operator-valued measure: (i) is the positivity con- 
dition, (ii) is a-additivity, (iii) is normalization. The instruments represent 
measurement procedures and their interpretation is as follows. íl is the set of 
all possible outcomes of the measurement ((Q, E) is called the value space) 
and the probability of obtaining the result uo G B {B G E), when before 
the measurement the system is in a state g (g G T(TC), g > 0, Tr{£>} = 1), 
is given by P(B\g) := Tr{2(B)g}. Moreover, let us consider a sequence of 
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measurements represented by the instruments Xi, X 2 , X n and performed 
in the natural order (X 2 after T\ and so on). We assume any time specifica- 
tion to be included in the definition of the instruments (Heisenberg picture). 
Then, the joint probability of the sequence of results ui\ G B\, uj 2 G B 2 , 
u n G B n , when the premeasurement state is g, is given by 

P(B 1 ,B 2 ,...,B n \g) = Tr{l n (B n )ol n _ 1 (B n _ 1 )o...ol 1 (B 1 )g} . (0.1) 

If we consider the conditional probability of the results lü 2 G B 2 , ui n G B n 
given the first result u% G Bi, we can write 

P(R R\R-n\ ~ P (Bi,B 2 ,...,B n \g) 
P(B 2 ,...,B n \B h g) = p^-j 

= P(B 2 ,...,B n \g(B 1 )) = Ti{l n (B n )o...ol 2 (B 2 )g(B 1 )} , (0.2) 

where we have introduced the statistical operator g{B\) representing the 
state after the first measurement, conditioned upon the result íü\ G B\. For 
a genèric instrument X and set B, the conditioned state g(B) is defined by 

!(í) -ïjïH = pS' (0 ' 3) 

Let us note that joint probabilities (|0.1|) preserve mixtures, by the linearity 
of the instruments: for g and a statistical operators and < À < 1, we have 

\P(B U ...,B n \g) + (l- X)P(B h ...,B n \a) = P(B 1 , . . . , B n \Xg + (1 - X)a) . 

(0.4) 

However, this property is not shared by conditional probabilities (10. 2 j) . by the 
definition itself of conditioning, and, therefore, the expression (10.3)1 for the 
conditioned state is not linear in the premeasurement state g, unless B = Q. 

Formula ()0.3)) can be interpreted by saying that we perform some mea- 
surement on a statistical ensemble of systems and select those systems for 
which the result u G B has been found. Then, (j0.3p is the state after the 
measurement of the systems selected in this way and depends not only on 
the result w6 5, but also on the perturbations due to the concrete measur- 
ing procedure and to the dynamics. If we perform the measurement, but no 
selection, we obtain by (iii) g(fí) = I(Q)g. By the definition of instrument, 
this quantity is linear in g and it is a statistical operator if g is a state. We 
can call g(fl) the a priori state: if we know the premeasurement state g 
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and the measurement X, g(Q) is the state we can "a priori" attribute to our 
systems, before knowing the result of the measurement. 

Let us consider now the case of the most fine selection when in ()0.3|) 
the set B shrinks to an "infmitesimally small" set du around the value u. 
According to the discussion above, the quantity 

X(dw)g 

= Tr{X(du)g} (0 ' 5) 

represents the state conditioned upon the result u G du. The quantity g{ui) 
is the state one can attribute to those systems for which the result u has 
actually been found in the measurement and for this reason we call it a 
posteriori state 

More precisely, a family of statistical operators {g(u), uj G Í2} is said to 
be a family of a posteriori states [3Ü] , for an initial state g and an instrument 
Z with value space £), if (a) the function u —>■ g(uj) is strongly measurable 
with respect to the probability measure 

f i g (B):= r It{X(B)g} = P(B\g) (0.6) 

for the observable associated with the instrument X and (b) VI" G B(7í), 
Vfi G S, 

/ Tr{Yg(u)}fx g (dLü) = Tr{F X(B)g} . (0.7) 

Let us note that by defmition the link between a priori and a posteriori states 
is given by 

g(Q)=X(Q)g= í g(u) ^(du) . (0.8) 

Let us stress that (|0.7|) defines the a posteriori states once the instrument 
X and the premeasurement state g are given. On the contrary, if g(uj) and \i Q 
are given for any g, (j0.7)l allows to reconstruct the instrument X. We shall 
make use of this in the following sections. 

Finally, let us note that there is no reason for g(u) to be a pure state if g 
is pure: it depends on the concrete measuring procedure. Roughly speaking 
£>(u;) is pure if one has some property of minimal disturbance, some ability of 
the measurement to give a maximum of information, we shall see various 
examples (g(u) pure and not pure) in the case of continuous measurements. 
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1 Counting processes 

The first class of continuous measurements which has been introduced in 
quantum mechanics is that of counting processes f2\ H] P Hü} IT71 HH} l2~%l IÏÏT] 

[3*lj . One or more counters act continuously on the system and register the 
times of arrival of photons or other kinds of partides. 

Let us consider the case of d counters. They differ by their localization 
and/or by the type of partides to which they are sensible and/or by their 
operating way... We can describe this counting process by giving the so called 
exclusive probability densities (EPDs) [7J |Ü]. The quantity Pí(0\g) is the 
probability of having no count in the time interval (to, t], when the system is 
prepared in the state g at time t . The quantity p\ Q (ji, t\, 32, t 2 ; ■ ■ ■ ; j m , t m \g), 
jk = 1, • • • , d, to < t\ < Í2 < • • • < t m < t, is the multi-time probability 
density of having a count of type j\ at time íi, a count of type & t time Í2, 
. . ., and no other count in the rest of the interval (t ,t]. Davies jlj (see also 
0-0) has shown that these EPDs can be consistently described in quantum 
mechanics in the following way. 

Let Coit) be a Liouvillian (the generator of a completely positive dynamics 
|2Ü] on T(TC)) and Jj{t), j = 1, . . . , d, be completely positive maps on T(7í). 
Let us introduce the positive operators Rj(t) on TC by 

R^-.^jjityi. (1.1) 

For any operation A on T(7í) its adjoint A' on B(TC) is defined by 

Tr{XAg} = Tr{gA'X}, Vg G T(H) , \/X e B(H) . (1.2) 



Finally, let S(t,t ), t > to, be the family of completely positive maps on 
T(Tl) defined by the equations 

|-S(Mo) = A(t)S(t,t Q ) , 5(í ,ío) = Id, (1.3) 
at 

A(t)g = £ (t)g-l-Ít{ R A t )iQ}- (1-4) 

Here {a, 6} = ab + ba and Id is the identity map on T(TC). Then, the 
quantities 

PÏ o (0\g)=Tr{S(t,t o )g}, (1.5) 
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PtoÜí» h',32, *a; ■ ■ ■ ; j m , t m \g) = Tr{S(í, t m )J jm (t m ) 

XS^tm-^Jj^itm-x) ■ ■ ■S{t 2 ,t 1 )J jl (ti)S{t 1 ,t )g} (1.6) 

(where ío < ti < t% < ■ ■ ■ < t m < t, jk = 1, • • • , d) are a consistent family of 
EPDs. 

The whole statistics of the counts can be reconstructed from the EPDs. 
For instance, the probability of m counts of type j in the time interval (í , íi], 
n counts of type i in the interval (t\, £2] is given by 

P(m,j, (t , íi]; n, «, (ti,t 2 )\g) = 

rt2 rr n n'2 rti rs m 

= dr n dr n _i • • • / dri / s m / ds m _i ■ ■ ■ 




dsi P%(j,si;j,s 2 ; . . .;j,s m ;i,ri;i,r 2 ; . . .;i,r n \g) . (1.7) 



In a similar way all more complicated joint probabilities can be constructed. 

One of the most significant problems treated by this theory is that of the 
electron shelving effect or quantum jumps. An atom with a peculiar level 
configuration and suitably stimulated by làser light emits a pulsed fluores- 
cence light with random bright and dark periods. It is possible to use Co(t) 
for describing the free atom and the driving term due to the làser and to 
use the operators Jj{t) for describing the emission process. Then, the full 
statistics of the fluorescence light can be computed and, in particular, the 
mean duration of the dark periods [HI-JI!!]- Other applications to quantum 
òptics of the counting theory described here are given in [35|-[37|. 

Let us now consider the problem of the a posteriori states. Our counting 
process can be considered as a stochastic process whose associated probability 
measure (uniquely determined by and (jl.fij) ) is concentrated on step 

functions. Let us consider t = as initial time. A typical trajectory u t up 
to time t is specified by giving the sequence (ji, í x ; j 2 , t 2 ; . . . ;j n ,t n ) of types 
of counts and instants of counts íi < t 2 < ■ ■ ■ < t n up to time t. Let uo t be the 
trajectory we have registered up to time t. Then, the conditional probability 
P(0, (t,t+t\\uJt', g) of no count in the interval (t,t+t\, given the state g at 
time zero and the trajectory ut, is given by 

p ( o, (t, í+ï]k g) = P j^' 1 ; t ! ; ··; ; / w ; t ;'f m 

P0\J1) r li • • • ) Jn, i>n\Q) 
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(cf. dü3)). By (fOjl . (H3D, dUSD, we obtain immediately that the probability 
(|1.8j) can be rewritten as 



P(0,(t,t+t\\iü t ;g) = Ti{S(t+t,t)g(t)} = P t + t (0\g(t)), (1.9) 

g{t) ^^S(t,t n )j jn (t n )S(t n ,t n ^)---j jl (t 1 )S(t 1 ,o)g, (í.io) 

where C(t) is the normalization factor determined by Tr{g(t)} = 1 (cf. ()0.2j) . 
(10. 3|) and (j0.5|0 . Similar results hold for the other EPDs conditioned upon 
some trajectory up to time t. Therefore, all conditional probabilities can be 
computed by (jl.5j) and (jl.fij) if one uses as initial state the expression (jl.lOJl . 
Equation (jl.K)j) gives the state of the system at time t conditioned upon the 
trajectory uj t up to time t (the a posteriori state). 

The interpretation of (|1.10|) is that, when no count is registered, the 
system evolution is given by S(t,to) and that the action of the counter on 
the system, at the time t in which a count of type j is registered, is described 
by the map <Jj{t). However, 5(í,í ) and J'jif) do not preserve normalization 
and the normalization factor C(t) is needed. This is due to the fact that they 
are the probabilities (jl.5|) and (jl.6j) which have to be correctly normalized 
and this is guaranteed by equations (|1.1|) and (jl.4|) . connecting Jj{t) with 
S(t, t ). According to this interpretation of ()1.10|) . the state of the system in 
between two counts is 

«M- Tf^'wMl - (L11) 
Tr{S{t,t r )g{t r )} 

where t r is the time of the last count and g(t r ) the state just after this count. 
If we denote Tr{Rj(t)g(t)} by (Rj(t)) t and differentiate (jl.lljl . we obtain 

^ = C (t)g(t) - \ £ {R 3 (t) - (R 3 (t)) tl g(t)} . (1.12) 

Moreover, if at time t r we have a count of type j, the state of the system 
soon after is 

,,,,,, JjiQgjtr) Jj(tr)Q(t r ) 

9{t r +dt) - Timtr)g{tr)} ~ {R . {tr))tr ■ (1.13) 

Now the typical trajectory Nj(t) (number of counts of type j up to time 
t), j = 1, . . . , d, of our stochastic process is a step function such that Nj(t) 
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increases by 1 if soon after time t there is a count of type j, otherwise Nj(t) 
is constant. Therefore, the Itò differential 

dNj(t) = Nj(t+dt) - Nj(t) (1.14) 

is equal to one if at time t there is a count of type j and to zero otherwise. 
This ffives (dJVj-(í)) = dNj(t). Moreover, the probability of more than one 
count in an interval dt vanishes more rapidly than dt, i.e. between dNj(t) 
and diVj(í), i ^ j, at least one of the two must be zero. Moreover, dNj(t) dt 
is of higher order than dt and has to be taken vanishing. Summarizing, we 
have the Itò table 

dNj(t) dJVi(í) = % dNj(t) , dNj(t) dt = 0. (1.15) 

By using these results we can rewrite (11.12)) and (|1.13|) in the form of a single 
stochastic differential equation in Itò sense ( dg{t) = g{t + dt) — g(t) ): 

dg(t) = C(Mt)dt + y t (|f|f - ( dN ^ - ( R Mtàt) , (1.16) 

C(t)g = C (t)g + é (Jj{t)g - \ {^(t), g}) , (1.17) 

{R,{t)) t = TiiR^git)} = TiiJ^git)} . (1.18) 

Indeed, when all dNj(t) vanish, (|l.lfj|l reduces to (jl.12)) : when one of the 
dNj(t) is equal to one all the other terms in the r.h.s. of (jl.lfij) are negligible 
and we obtain ()1.13J) . Equation ()1.16|) was firstly obtained by quantum 
stochastic calculus methods in |23 123 EE EH] ■ 

Formula ()1.16j) is the equation for the a posteriori states in the case of a 
counting measurement: it determines the state at time t depending on the 
(stochastic) trajectory up to time t. Let us stress that we know the solution 
of this equation: it is the state (jl.lOj) . In any case, it is very useful to have 
the differential stochastic equation f)l. 16|) as we shall see in the rest of this 
section and in section El 

Let (dNj(t))(üü t ) be the mean number of counts of type j in the interval 
(t,t+dt\ conditioned upon the trajectory u t up to time t. Because probabil- 
ities of more that one count in a small interval are negligible, we have 

(dNj(t))(u t ) ~p t j dt (j,t\g(t))dt~Tr{J J (t)g(t)}dt= (R^dt. (1.19) 
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In other words the quantities (Rj(t)) t dt appearing in are the a poste- 

riori mean vàlues of dNj(t). Moreover, the differentials 

dMj(t) = dNjit) - (Rj(t)} t dt, (1.20) 

appearing in (jl.lfíj) . together with the initial condition Mj(0) = 0, define the 
a posteriori centered processes Mj(t), called innovating martingales. 

Equation (|1.16|) is non-linear, but it is mathematically equivalent to a 
linear one. Let us introduce an arbitrary stochastic real factor c(t) and define 
the trace-class operator <p(t) := c(t)g(t). If we know <p(t) we can reobtain 
the state g(t) simply by normalization. The factor c(t) can be chosen in such 
a way that (p(t) obeys a linear stochastic differential equation; moreover, 
this choice is not unique. We shall do this in a very convenient way: a new 
linear stochastic equation is obtained giving both the a posteriori state g(t) 
and the EPDs (jl.5|) and (cf. jSH])- Let ip(t) be a trace-class operator 
depending on the trajectory oj t and defined by cp(t) = S(t,t r )(p(t r ) if t r is 
the time of the last count and by <p(t r +dt) = Tjj{t r )(p(t r ) if at time t r there 
is a count of type j; r is an arbitrary parameter with dimensions of a time, 
which disappears from the physical quantities. For initial condition we take 

<p(0) = Q. 

By the definition of (p(t), the quantity 

c(í)=TrMí)} (1.21) 

gives the EPDs (jl.5|) and in the case of a trajectory uj t containing no 

jump we have 

c(t) = P í (0|e) (1.22) 

and in the case of a trajectory with a jump of type j\ at time íi, of type 
j m at time t m we have 

c(t)=r m pl(j 1 ,t 1 ;...;j m ,t m \g). (1.23) 

Moreover, by the definition of <p(t), c(t) and g(t) the a posteriori state is 

g(t) = <p(t)/c(t) . (1.24) 

In the same way as for g(t), we can obtain the stochastic differential 
equation for f(t), which turns out to be (cf. [2E1, equation (??20)) 

dip{t) = [Ca(tMt) - \ ER-(í), ?(t)}]dt + YfrJ&Mt) - vmN^t) . 

Z 3=1 3=1 

(1.25) 
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By using Itò's calculus for counting processes it is possible to verify triat 
indeed (TT2TJ) . (Q3J) and (fH?ïïJ) are equivalent to (fTTïï|) . Equation (fOïïj) 
determines all the probabilities via f)1.21[) — (jl.23j) and the a posteriori states 
via (jl.21j) and (jl.24J) . Equation ()1.25[) is linear, once a realization Nj(t), 
j — 1, • • - , d, of the process is given. However, let us note that the statistics of 
Nj(t) depends in its turn on the premeasurement state g, as shown by (jl.5|) . 
(jl.òj) . (|1.19|) . The possibility of finding a linear equation mathematically 
equivalent to means that g(t) is linear in g up to a normalization 

factor, as suggested by (|0.5|) . 

Let us stress that in general equation (jl.l6j) does not transform pure 
states into pure states. This simply means that in the course of time we 
loose information due to some dissipation mechanism, for instance the system 
interacts also with some external bath or similar things. In any case, the 
situation in which pure states are preserved is particularly interesting. This 
is the case |H1 121 123 I2Ü when 

C (t)g = -i[H(t), g] , Jj(t) = ZjWgZjity , (1.26) 

where Zj(t) and H(t) are operators on 7í, H{t)^ = H(t). Then, we have 
Rj(t) = ZjltyZjit) and 

C(t)g = -i[H(t), g] + E (z^gZjitY - \ {z^Z^t), g}^j . (1.27) 

Then, f|1.16|) becomes 

dg(t) = -i[H(t), g(t)]dt - \ è {Z 3 (tyZ 3 (t) - (Z^YZ^t, g(t)}dt + 



By using Itó formula ()1.15|) . one cane prové that g(t + dt) 2 = g(t + dt), if 
g(t) 2 = g(t); therefore, (jl.28j) transforms pure states into pure states and it 
is equivalent to a stochastic differential equation for a wave function. Indeed, 
let %fj{t) e TC satisfy the "a posteriori Schrodinger equation" (2H1EH1 

dm = [-iff(f)-~é(z,(f)tz^ 

i=i 

+ E[ . Zj{t) =-i|#)djyj(f). (i.29) 
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with (Z j (í) t Z j (t)) í = (jl){t)\Zj{tyZj{t) ip(t)); then, by O) one obtains that 
g (t) = \ip(t)){ijj(t)\ satisfies (|ÒHJ). 

It is interesting to note that in between two counts (when dNj(t) = 0) 
(I1.29J1 becomes a nonlinear Schròdinger equation of the type studied, for in- 
stance, in 001 El- However, this equation has a quite different interpretation 
in the quoted references, where the problem is to find evolution equations 
compatible with the Hilbert space structure and preserving "properties" in 
the sense of quantum lògic. 

Now we have the a posteriori states defined by (J1.16J) and a probability 
measure on the trajectory space, which is implicitly defined by the EPDs 
(II. 5j) and Therefore, we can reconstruct the instruments associated to 

our measurement by means of (|U.7|) . As in j2Hl EH] , we shall do this by using 
the notion of characteristic operator, a concept introduced in [TÜJ [TÜJ , and 
Itò formula for counting processes. 

Let / be any function of the trajectòries of our stochastic process and let 
us denote by (/) st the mean value of / with respect to the measure associated 
to the EPDs (HISD and (fO)|l . The quantity 

= (exp{i£^%( S )dA^( S )}^ (1.30) 

is called the characteristic functional of the process. Here k(s) is a test 
function, i.e. kj(s) is a real compact support C^-function on (0, +oo). $t[&] 
determines uniquely the whole counting process up to time t: roughly speak- 
ing $t[&] is the Fourier transform of the probability measure of the process. 
More explicitly ^7], we have 

_ 00 00 rt ft m 

$ í [fc]=P t (0| í? )+E E dt m dí ro _i..- 

■■■j^ dt 1 exp^iJ2kjX t i)jPoÜi^iij2,t2]...-J m ,t m \g). (1.31) 
Let us set now 

V t [k] =exp{iX: [k^dNjis)}. (1.32) 
j=i J0 J 

According to í)0.7|l . we can write 

(V t [k]g(t)) st = g t [k]g, (1.33) 
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— * 

where g(t) is the a posteriori state at time t and Gt[k] is an operator on T(H) 
which represents the "functional Fourier transform" of the instrument X t 
associated to our measurement up to time t. The quantity Gt[k] can be called 
characteristic operator and it is the operator analogue of the characteristic 
functional of a stochastic process JI] JB]. By the normalization of g(t) and 
(Pül) and (fPïïl) . we obtain 



$ t [k]=Tr{g t [k] e }. 



;i.34) 



An equation for Gt[k] can be found by differentiating f)1.33|) . The differ- 
ential of g(t) is given by (jl.lfíjl . while the differential of V t [k] is 



dV t [k] = V t [k]\j2( 



1 dNAt) 



■3=1 



;i.35) 



This formula can be easily obtained from (jl.32|) by expanding the exponential 
and using (jl.lòj) . By using the formula 

d(V t [k]g(t)) = (dV t [k]) g(t) + V t [k](dg(t)) + (dV t [k)) (dg(t)) , (1.36) 

where the Itó correction (dV t [k\j (dg(t)) has to be computed by means of the 
Itò table f)1.15|) . we obtain 

d(V t [k]g(t)) = \/ í [fc]{/:(í) í? (í)dí + X:(e i ^W-l)^(í) í? (í)dí + 



d 

+ E 

i=i 



(dN 3 (t) - (R^tdtÚ^) 



Now let us take the stochastic mean of (jl.37|) . We compute this mean 
in the following way First we take the mean with respect to the probabil- 
ity measure on the future (with respect to t) conditioned upon the given 
trajectory. All the quantities in the r.h.s. of (jl.37J) depend only on the past 
(they are adapted), but the quantity dNj(t), whose a posteriori mean value is 
just {Rj(t)) t dt (equation ()1.19|) ). Therefore, the last term in (jl.37j) vanishes. 
Then, we take the mean value also on the past and, by ()1.33|) . we obtain 



j t Gt[k] = JC t (k(t))Ç t [k] 



;i.38) 
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)C t (k(t)) = C(t) + £ (e^ {t) - l) J 3 {t) ■ (1-39) 
Together with the initial condition 

G$\ = Id (1.40) 

(which follows from the definition (jl.33p ). equation (|1.38|) determines uniquely 
Gt[k] and implicitly the instruments on the trajectory space. This kind of 
equations has been obtained for the first time in [T3] . 

If no selection is made according to the results of the measurement (let 
us say: the results are not read), the state of the system at time t will be (cf. 

(EBD) 

<t(í) = <*?(*)>*; (1-41) 

a(t) is the a priori state for the case of the continuous measurement described 
in this section. According to (021) . IfDEty . (OBI) . (fOUjl . we have that the 
a priori states satisfy the quantum master equation 

±a(t) = C(t)a(t) , (1.42) 

with the new Liouvillian ()1.17j) : the unperturbed Liouvillian Coit) corrected 

by the measurement effect term Yfj=i {jj{t)Q ~ \ {Rj(t)i í?})- The fact that 
we have obtained a linear equation for the a priori states is due to linearity 
and normalization of the instruments (cf. ()0.8|0 . 



2 An example of counting process: 
a two— level atom 

Let us consider an example of counting measurement on the simplest quan- 
tum system: a two-state system, described in the Hilbert space TL = C 2 . 
We can think of a two-level atom, an unstable particle, a spin... While the 
general case could be handled, for concreteness we treat a two-level atom 
with pumping and damping. This section has to be considered just as an 
illustration of the theory developed before. 

The (time independent) unperturbed Liouvillian is given by 

C g = -^uj[a 3 ,g] + JoQ-^{Ro,q} , (2.1) 
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JoQ = cr+Qcr- + A_ (T-Q(J + , (2.2) 

#o = JÓ<ro = \ A+(o - + ~ A-((To + 0-3) • (2.3) 

Here u > 0, A± > 0, Oj, i = 1,2,3, are the Pauli matrices, Oo is the 2x2 
identity matrix and a± = ~(oi ± io 2 ). 

We consider a single counter (d = 1) and take 

J 1 g = \ 1 (T_Qa + , Xi>0; (2.4) 

the map J\ describes the emission of photons (or other types of partides), 
which are then counted by some electrònic device. In the present case, the 
rate operator (|1.1|) is 

Ri = J(cro = Ai o + o_ = -Ai (0-0 + 0-3) (2.5) 
and the generator C ()1.17j) of the full dynamics is 

Cq = - 1 -uj[<t 3 , q] + £ (jjQ - \ {Rj, £?}) ■ (2.6) 

We can intèrpret the terms with A + as pumping, the terms with A_ as 
incoherent damping and the terms with Ai as electromagnètic decay; T = 
Ai is the electromagnètic transition rate. If A + = 0, we can intèrpret the 
system as a Wigner atom (or another unstable particle). In this case the 
electromagnètic transition rate is T — A_ + Ai; A_ 7^ means that not all the 
photons are collected by the photocounter; e — X\/ (A_ + Ai) is the efficiency 
of the counter [T7] . 

In order to perform computations, it is convenient to represent selfadjoint 
trace-class operators (p as 

p=Í(ar + C<T+ + CV_+foO, c,eelR, (GC. (2.7) 

£ + |C| ) an d it is a density matrix if 

also c = 1. 

Let us consider (jl.25j) and represent (p(t) in the form (|2.7|) with c — * c(t), 
C — *■ C(^)> £ The stochastic equation (|1.25p . choosing r = A x , 

becomes 

dc(t) + í Ai[c(í) + i(t)W = \ W) - c(t)]dN(t) , (2.8) 
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d£(t) + 



2«-7ïAi)e(í) + (a + -Ai)c(í) 



dí = -±[c(í) + 3£(í)]diV(í) 



(2.9) 

dC(t) + (io; + «)C(*) dí = -C(í) dN(t) , (2.10) 

where k = ^(A + + À_ + Ai), a = À_ — À+. It is convenient to rewrite f)2.8j) 
and ()2.9|) in terms of the stochastic parameters 



7T (í) = i (c(t) - eco) , = i (c(o + ew) ; (2.H) 



this gives 



dTTo(í) + [^ T 7To(í) - «iTTiíOldí = K (í) - 7T (í)]diV(í) , 

dTTi(í) + [pt^iit) - « T 7T (t)]dt = -7n(í)dJV(í), (2.12) 

where ^ — k-j- — A+, K| = À_, //i = À x + A_. 
The solution of (|2.1üj) is very simple: 

where t x is the instant of the first jump of N(t). About (|2.12|) . let us denote 
by 7Tj(t\a,b) the solution of (|2.12|) with dN(t) = and initial conditions 
7r (0) = a, 7Ti(0) = b. Then, the solution of the stochastic system ()2.12j) is 

A) \ ^(í — í r |7ri(í r ), 0), if í r < t < t r+1 , r > 1, K > 

where t r are the instants of the jumps of N(t). 

By í)1.21|l . (|1.24j) . ()2.7|) . the matrix elements of the a posteriori state g(t) 
are given by 

(1|é?(í)|1> = Tr{~(a + CT3)Q(t)}=Tr 1 (t)/c(t), 

(0\g(t)\0) = Ti{^(a -a 3 )g(t)} = n (t)/c(t), 

(l\g(t)\0) = Tr{a^(í)} = C(í)/[2c(í)], 

(0|e(í)|l> = Tr{a +6 (t)} = ((ty/[2c(t)\, (2.15) 

with c(t) = 7To(í) + 7Ti(í). Equations (|2.13jl - (|2.15jl shows that at a jump of 
N(t) the system surely goes into the ground state, because ( = and iï\ = 0, 
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and that for t > t% the system is surely in a mixture of ground and excited 
states, because £ = 0. The EPDs are implicitly given by c(t) = ir (t) +7r 1 (í), 

r = Ar 1 , <B, (EU- 

Just as an example let us discuss the case of the Wigner atom (À+ = 0). 
Equations ()2.14j) become 



7n(í) 



' *o(0) + ^(l-e- 2 « í )7r 1 (0), ifí <íi, 
exp[-2KÍ 1 ]vr 1 (0), ifíi<í<í 2 , (2.16) 

k 0, ifí >í 2 , 



with k = |(A_ + Ai), 7T (0) + 7Ti(0) = 1. Equations (EHoT) give = |0)(0| 
for t>t\. after the first registered emission the atom is in the ground state. 
Finally the EPDs are 

P*(0\g) = tto(O) + -L (A_ + A ie - 2KÍ ) 7n(0) , (2.17) 



2k 

PoUiMq) = A i exp[-2«í 1 ]7n(0), (2.18) 
Pl(juti;...;j m ,t m \g) = 0, m > 2 . (2.19) 

These equations say that there is at most a count, as it must be because 
there is no pumping. 



3 Diffusion processes 



In the classical case Gaussian diffusion processes can be obtained from Pois- 
sonian counting ones by centering and scaling. Similarly, in the quantum 
case we can obtain some kind of "quantum diffusion measuring processes" 
from the quantum counting processes of section ^ 

Let us take the maps Jjit\ describing the action of the counters, of the 
following form: 



1 



j 3 (t) Q = z j (t) + -f j (t) g z s {ty + -f^ty 



i 



(3.1) 



where the Zj(t) are operators on TC, the fj are complex functions and e > is 
a parameter which we want to make vanishing at the end. Moreover, instead 
of £o(t) we take as unperturbed Liouvillian the expression 

Cl(t)g = C (t)g + J- é [ifAtyZAt) - W)Zj(tV, g] ■ (3.2) 
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Then, the generator C(t) of the a priori dynamics (cf. (|1.17|) and Q1.42]) ) 

becomes 

C(t)a = C (t)a + ^(ZjQaZ&y - \ {z^Z^t) , a}) . (3.3) 

The expression ()3.2|) has been assumed in order to have C(t) independent of 
the parameter e. Physically, the structure f)3. lj) — f)3.3j) is related to heterodyne 
detection 

Moreover, we make a linear transformation on the outputs: we call Y?(f) 
the new observed processes, related to the old processes Nj(t) by 



dK/(í): = ^dA^(í) - - l/^í)! 2 dí: 



(3.4) 



this means that we rescale the outputs and subtract a known deterministic 
signal. Then, by ()3.4j) and (jl.l5|) we obtain 



dY/(t) dYf(t) = e 2 b lj dN^t) =e5 i3 dF/(t) + 8ij \f 3 {t)\ 2 ^ • 



(3.5) 



In order to have the characteristic operator associated to this new pro- 
cesses, we have to rescale the test function k(s), appearing in (jl.30|) . (|1.34|) . 
(|1.38)l — (jl.40jl . by changing kj(t) into ekj(t) and we have to shift the mean 



vàlues of sNj(t) as in ()3.4|) by adding to )C t (k(t)) the term — ^ J2j kj(t) \fj(t)\ . 
The final result is that the generator JC t (k(t)) of the characteristic operator 
Gt[k] becomes 

K t {k{t)) e = c{t) Q + j^\-\k 3 {tf\fM 2 Q + ^\f^rz 3 {t)Q 



+ fi{t)QZj(tf\ 



ZAt)gZAt) 



1 

+ - 

e 



iekjit) [fjityZ^g + fj^gZ^ 



(3.6) 



Also equation (jl.lfij) for the a posteriori states can be expressed in terms 
of the new processes Yf(t). By (JTÏJ), Q| . (jOl . (I3~ÏÏ|) and (Q, we obtain 



df?(t) = C^g^dt + ^^Z^gZ^tY -eiZ^tyZ^tgit) 
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;e 2 <Z i (í)tZ i (í)> t + e/ i (í)·<Z i (í)>* + e/ i (í)(Z i (í)t> t + |/ i (í)| a ; 
d!7(í)-e(Z i (t)tZ i (í)) í dí-/ i (í)*(Z i (í)) í dí 



-^•(íK^t^dt , 

where, for any operator X on 7í, (X) t is defined by 

(X) t = Tr{Xg(t)}. 



(3.7) 



(3.8) 



Moreover, from (|3.4j) . ()1.19j) and ([3.1jh we have that the a posteriori mean 
vàlues of dY ? e (t) are given by 

(dY/(í)>(u, t ) = [/ J (í)(^(í)t) í + / J (í)*(^(t)) í ]dí + 5(^(í)t^·(í)) í dt. (3.9) 

We assume \fj(t)\ ^ 0, Ví. From IJ3.5j) — (|3.7j) and ()3.9|) . it is apparent that 
the limit e l exists. In this limit we obtain that the characteristic operator 
is given by (jl.38|) and (jl.4()|) with generator 

K t {k{t))g = C(t)g 

+£{4 Htf \fM 2 q + ikMfAtyzMe + /;(*) ^(t)t]}(3.io) 

3=1 

By setting Yj(t) = lim e j_o Yjif), the equation for the a posteriori states be- 
comes 

dg(t) = C(t)g(t)dt 

+ EifÀtriZjit) - (Z 3 {t)) t ]g{t) + / i (í)^(í)[Z,(í)t - (Z s (t)%)} 



3=1 



X 



dYJt) - fJtYiZ^tdt - fAt){ZAtf) t di\ . (3.11) 



Moreover, the a posteriori mean value of dYj(t) becomes 

(d^(í))(w t ) = 2Re[/ i (í)*(Z i (í)) t ] dt (3.12) 
and the processes Mj(t), defined by 

dMj(t) = dYj(t) - 2Re[f j (t)*(Z j (t)) t ] dt , M,(0) = , (3.13) 
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are again innovating martingales. Finally, the multiplication rule for the 
differentials dYj(t) is the limit of (|3.5j) under e { 0; by taking into account 
also the second of equations (jl.l5|) . we have the Itò table 

dY 3 (t) dY t (t) = Sjí \fj(t)\ 2 dt > dY j(t) dt = • ( 3 - 14 ) 

By the procedure we have followed, it turns out that also the connection 
between a posteriori states g(t) and characteristic operator Ç t [k] given by 
()1.33|) continues to hold, but now g(t) satisfies ()3.11|) . Gt[k] satisfies ()1.38|) . 
(jl.40j) with generator given by (j3.10|) and Vt[k] is given by 



(3.15) 



V t [k] = expli^t fk^àY^s) 
L j=1 Jo 



Alternatively, equation (jl.33j) can be proved by taking the stochastic differ- 
ential of both its sides, as done in the case of counting processes. 
By taking the mean value of ()3.12j) on the past, we obtain 

f t (Y 3 (t))st = T±{MY Z Àt) + fj(t) Z AtfW)} , (3-16) 

where a(t) are the a priori states satisfying equation (jl.42|) with Liouvillian 
(13. 3J) . The same result can be obtained by functional differentiation of the 
characteristic functional Tr Qj, T > t, with respect to kj(t) [11. 

Equations ()3.15|) and ()3.16|) show us two things. First, our continuous 
measurement gives the statistics of the generalized derivatives jlS] Vj(t) = 
Yj(t) (or of the increments dYj(t)) more than the statistics of the Yj(t) them- 
selves. The same was true in the case of counting processes, but in that 
case this difference was irrelevant, because we had the natural initial con- 
dition iVj(O) = 0. Second, ()3.16|) can be interpreted by saying that yj(t) is 
the output of a continuous measurement of the quantum observables (selfad- 
joint operators) fj(t)*Zj(t) + fj(t)Zj(tY , which are in general noncommuting 

[ÏÏ3-[I21E2H25!- 

Measuring processes defined by a characteristic operator with generator of 
the type ()3.1Ü|) were introduced in and equation (|3.11|) was obtained 

by quantum stochastic calculus methods in |23]-j2ZlEHlllll· By linear trans- 
formations on the outputs, the most general diffusive case can be reached; 
moreover, by taking prescription (|3.1|) only for a subset of the J~j, mixtures 
of diffusive and Poissonian contributions can be obtained (2E]- Itò equations 
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for the a posteriori states in the jH3 HüJ E7j purely diffusive case have been 
considered also by Diósi |%5]-|T7j. 

As in the case of counting processes there exists a (not unique) linear 
stochastic equation mathematically equivalent to ([3. II)) . For instance, let 
<f(t) be a trace class operator satisfying the equation 38J 

d 

d<p(t) = C(t)(p(t)dt + 

and set c(t) := Tr{(p(t)}. Then, by Itó's calculus one shows that g(t) = 
ip(t)/c(t) satisfies (|3.11jl . To the linear equation 1)3.17)1 the same comments 
apply as to (jl.25|) . 

In the case of an unperturbed Liouvillian of a purely Hamiltonian form, 

C (t)g = -i[H(t), e ], (3.18) 

equation ()3.11j) transforms pure states into pure ones; for proving this it 
is sufficient to show that g(t + dt) 2 = g(t + dt) if g{t) 2 = g(t). In this 
case, which we can call of complete measurement, ()3.11j) is equivalent to a 
stochastic differential equation for a wave function, as in the case of counting 
measurements. Indeed, let ip(t) 6 7i satisfy the "a posteriori Schròdinger 
equation" [23 Ej 

dm = -{iH(t)+litiz J (tyz J (t)-2(z J (ty) t z J (t) + \(z,(t)) t \ 2 ]}mdt 

j=i 

-miZ^tdt], (3.19) 

with (Zj(t)) t = (ip(t)\Zj(t) ip(t)); then, by Itò's calculus, one obtains that 
g(t) = \^{t))(i(>(t)\ satisfies (TTITIt . 

It is interesting to note that stochastic equations of the type of í)3.11j) and 
(J3.19j) . with fj(t) = 1, have been appeared in the literature also in connection 
with dynamical theories of wave-function reduction jH^jHS- The idea is 
that the wave-function reduction associated to a measurement is some kind 
of stochastic process and an equation of the type of (j3.19j) is postulated. 
Apart from the different interpretations, another important difference is that 



àYj(t) (3.17) 
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in the dynamical reduction theories the noise comes from outside, while for 
us it is determined by the system itself. 

Sometimes it is useful to have at disposal a complexified version of dif- 
fusion processes. Let us consider the case of an even d. By redefining d 
and the index j, the sum appearing in (j3.10|) and (|3.11j) can be reorganized 
as a double sum over A, A = 1,2, and j, j = l,...,d. Then, we take 
= 1, f 2j (t) = i, Zij(t) = Z 2j (t) = Zj (t) and set K j (t) = hjjt) + ik 2j (t), 
Wj(t) = \ ÍXij{t) + iY 2j (t)). Then, (Q, dSIïnD (EH2D , §M become 

C(t)a = C (t)a +2^fe(í)aZ i (í)t - \ {z^Z 3 {t) , *}) , (3.20) 

K t {K{t)) e = C(t)g + £{-1 \ Kj (t)\ 2 e + ^{tyZ^Q + ^(í) ^(í) f ]} , 

(3.21) 

df?(t) = £(í)^(í)dí + 2^{[Z i (í)-(Z i (í))^(í)[dW i (í)*-(Z i (í)t) t dí] 

i=i 

+ [d^-(í) - (Z i (í))tdí]p(í)[^·(í) t - (Z^)*]} , (3.22) 

dW,-(t) àWi{t) = , dW,(t)* dWi{t) = - 5 Si dt , dWj(t) dt = , (3.23) 

(d^-(í))(a; í ) = (Z 7 ·(í)) í dí. (3.24) 
By taking the mean value of (|3.24j) on the past, we obtain 

j t {W 3 {t)) st = Ti{Z 3 {t)a{t)} 1 (3.25) 

which allows to intèrpret the equations above as describing a continuous 
measurement of the noncommuting, nonselfadjoint operators Zj(t). Filtering 
equation ()3.22|) for linear systems (quantum oscillators) was introduced in 

[211 123- 

4 An example of diffusion process 

Let us close the paper with a simple example of the theory developed in 
SectionEl in the complexified version (J3.2(J)) - (J3.25|) . A real-valued Gaussian 
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example for an observed particle in a quadratic potential can be found in 
ES]- We consider a single-mode field in a cavity and with a source, 

H(t) =ua ] a + g{t)a ] + g{t)*a, uj > , (4.1) 

interacting with a thermal bath, 

C (t)g = —i[H(t), g] + Xi([ag, a f ] + [a, gaJ]) + X^([aJg, a] + [a 1 ", ga]) , (4.2) 

Àj, À| > 0, and subjected to the measurement of a single complex observable 
(d = 1) proportional to the annihilation operator, 

Z = r?a, ?7<EC. (4.3) 

The fact that Z is proportional to a means that we are considering a passive, 
purely absorbing detector. 

By scaling the output in such a way that we have exactly a measurement 
of a (dW(t)/r] -> dW(t), rj*K(t) -> k(í)), equations (|ÏÏ^Üj) - (l3~23j) become 



£(í)cr = C (t)<j + |??| 2 (jacr, a 1 ] + [a, era 1 ' 



i 2 r + 

/C t (/í(í) , K{t))g = £(t)g — - \K[t)/r)\ g + i K[t)*ag + K{t)ga/ 



dg(t) = C(t)g(t)dt + 2\r]\ 2 {[a- (a) t ]g(t) [dW(t)* - (a)) t dt 
+ [dW(t)-(a) t dt] g(t) [at-(at) t ]}, 



(4.4) 
(4.5) 



(dW(t)) 2 = 



d(W(t))(u t ) = (a) t dt. 



dW{t) dt = 
■(W(t)) 8t = Tr{aa(t)}. 



(4.6) 
(4.7) 
(4.8) 



|dW(í)| = dí/(2|?7| 
d_ 
dí 

Equations f|1.38j) . with generator (|4.5|) . and ()4.6|) can be solved by anti- 
normal ordering expansion of trace-class operators. Let us define on T(7í) 
a "tilde" operation by tp e T(7í) — > 



^(r,0=Tr{e- i «*>e- i « at } 
which can be inverted by 



<P = ~ I d 2 ^é^ a e^(p(e,0 

71 



(4.9) 



(4.10) 
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Let us set 

<p(t) = (p(n*, «; t) = Q t [K*, k] q ; 
then (jl.38j) and ()4.5|) give, by Standard computat ions, 



(4.11) 



ílü + - r cd* + (h}-~r)çd- K{t)*o* - K (t)d 



2^\t\-9{t)e + g(tyt--\K(t)/7,\*\<p(e,£,t), (4.12) 



where d = d/d£, d* = d/dC an d 



r = 2 



+ X l - Af) . 



(4.13) 



We suppose T be strictly positive. 

If the initial condition is "Gaussian" . 



g(t,0 = exp 



(4.14) 



then <y9 maintains this structure at any time. Indeed, by inserting 



ï(e , e; í) = ex P {-i [«i) + - - [t 2 <*) + m*)*] - i^i 2 m - *(*)} 

(4.15) 

into ()4.12j) . we obtain the differential equations for the coeficients (/ is real): 

b(t) = 

c(t) = 

d(t) = 
f(t) = 

h(t) = -m(t)*b(t) -m(t)c(t)* + - \K{t)/r]\ 2 . (4.16) 



+ - rj 6(í) + ÏK{t)*d{t) + i«(í)/(í) - ip(t) . 

ílü + i r) c (í) - ík(í)M(í) - í«(í)/(í) - uj(í) 

-(2k; + r)tZ(t) , 
-r/(í) + 2A T , 

-ík(í)*6(í) - ík(í)c(í)* + - \K{t)/r]\ 2 . 



The solution of these equations can be easily written down. 

The characteristic functional of our generalized process [13] W(t) is given 
by (see OI) and (OÜ1) ) 



k] = Tr{^(/í*, k; í)} = £(0, 0; í) = exp[fc(t)] 



(4.17) 
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with 



h(t) = -i / ds [«(s)*a(s) + n(s)a(s)*} + [ dsds' \k(s)*k(s')Ai(s, s' 
Jo Jo 1 

+ i «(s)/c(s')A 2 (a, s')* + \ «(«)*«(«')*A 2 (s, s') 



(t) = e-^+s^ao - i Pdsgfàe 

Jo 



(iw+r/2)(t-s) 



(4.18) 
(4.19) 



r5(s - s') + t?(s - s') e 



-(iíj+r/2)(s-s') 



-H?(/-s) e (i ^ r/2)(s '- s) C(s), 
C (*) = 2 -^+U- 2 -p) e-* 



A 2 (s, s') = e 



-(iw+r/2)(s+s') 



/V 



(4.20) 

(4.21) 
(4.22) 



where $ is the usual step function. <&t[k*, is the characteristic functional 
of a Gaussian complex process with covariance ()4.2()j) . (|4.22jl and (a priori) 
mean vàlues 



K=K*=0 



Tr{aa(t)} = a(t) 



(4.23) 



The a priori states are given by a(t) = Gt[0) Q or <x(£*, £; í) 



, • B y (11351) and iftTiïijl we obtain 



K=K*=0 



<?(£*, £; t) = exp { -i iCa(t) + ce] - - C 2 e 
This gives 

Tr {«tfldfí)} = C(t) , Tr{a 2 (j^j| = exp 



r*2 -(2iu)+r)t 



-2[L·o + -Y)l 



(4.24) 



A*o- (4-25) 



Note the links between the covariance (J4.25j) of the a priori states a(t) and 
the covariance (jOüjl . (J4~22|) of the process W^í). 
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By the "tilde" transformation (|4.9|) . we can solve also the equation for 
the a posteriori states (|4.6|) . From (|4.fí|) . (|4.9|l . (|4.1U|) we obtain 



^liu + -r)ed* + (iu--r)td-2\ í \tf-g(t)e 



+ 9(tn] ~g(e, C; t) dt + 2 |r/| 2 { [dW(í)* - (a^tdí] - (a) t ] 
+ [dW(t) - (a) t dt] [id - ( a t) t j }e(r, C; í) • (4.26) 

This equation can be rewritten in terms of the stochastic function 

Z(t, £;*): = -ln (4.27) 

By using Itó's formula dg/g = —dl + |(d/) 2 , which in turn implies (dg/g) 2 = 
(dl) 2 , and Itó's table (l4~7j) . we obtain 



dl 



2 \rj\ 2 d*l dl - [iu + - T Cd*l + iw - - T ) Çdl 



+ 2A T K| 2 + g(t)C - - 2 |rç| 2 |(a) t | 2 ]dí 

+ 2 |r/| 2 \(id*l + (a) t ) dW(ty + (idl + (a%) dW(t) 



(4.28) 



With the initial condition (|4.14|) the solution of ()4.28|) remains quadratic 
in £ and £*. Indeed, let us write 

kc, 6 *) = i[r (o> t + £(at) t ] + i + em*\ + iei 2 , (4.29) 

where u(t) > 0; the term independent of £ is lacking because of normalization 
of g(t) and the linear term must have just the form we have written because 
(a)t is the a posteriori mean value of W(t). By inserting ()4.29|) into (|4.28J) 
and equating the coefficients of the same order in £ and £*, we obtain 



d(a) t + 



icü + -T)(a) t + ig(t) 



dt 



2 \ V \ 2 {n(t) [dW(t)* - (a + ) t dí] + z/(t) [dW(t) - (a) t dt]} , (4.30) 
d 



dt 



fi(t) + (2iu + T) fi(t) = -4 \t]\ z n{t)v(t) 



±u(t) + Tv(t) = -2 M 2 (|/i(t)| 2 + v{t) 2 ) + 2A T , 



(4.31) 
(4.32) 
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with (aJ) t = (a)í and the initial conditions (a) Q = a , /i(0) = /x , z/(0) = u . 

In the = 0, we obtain =0 (the stationary solution of (|4.3ip ) 

and (|4.32|) becomes 

^u(t) + Tu(t) = -2 \r]\ 2 v{tf + 2A T , (4.33) 
which is Riccati's equation and has the stationary positive solution 



r 



4 



1 + 16 \r)\ 



r 2 ) 



1/2 



(4.34) 



Equations ()4.3()j) (for /x = 0) and (J4.33j) were obtained for the first time in 
|2*H |23 121] as optimal filtering equations for linear systems. 

After a transient any memory of the initial condition is lost. The charac- 
teristic functional is given by (j4.17j) and (j4.18j) with a priori mean vàlues 

rt 



a(t) = 

and covariance A 2 (s, s') = 

1 



A 1 (s,s') 



-(io;+r/2)(í-s) 



2A 



g(s) ds 



5( s - s >) + ^l e -(m\s-s'\ e -Hs-s>)_ 



(4.35) 



(4.36) 



The a priori states are given by 



érooCr , £; í) = exp { -i [f a(í) + - ^ |£| 2 



while the a posteriori states are 

f?oo(C, 6 *) = ex P {-i [r(o)í + £(a f )í 
with a posteriori mean vàlues 



Kl 2 } 



( Q )< 



exp 



iu + - T + 2 |r/| z/oo ) (í - s) 



-hj(s) ds + 2|?7| 2 z/ 0O dW(s) 



(4.37) 



(4.38) 



(4.39) 



Note that > 2X^/T for Aj > and = for A-f = 0. In this last case 
the asymptotic a priori and a posteriori mean vàlues coincides (at(t) = {a) t ) 
and the same holds for a priori and a posteriori states 

a 00 (t) = g 00 (t) = \a(t))(a(t)\ (4.40) 

where \a) denotes the usual coherent states and a(t) is given by ()4.35|) . 
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